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We derive a simple analytical expression for the level correlation function of an integrable system. 
It accounts for both the lack of correlations at smaller energy scales and for global rigidity (level 
number conservation) at larger scales. We apply our results to a rectangle with incommensurate sides 
and show that they are in excellent agreement with the limiting cases established in the semiclassical 
theory of level rigidity. 

I. INTRODUCTION 

In the quantum limit, classically integrable systems are often described as those with "level bunching" or, alterna- 
tively, absence of level correlations, while classically chaotic ergodic systems are described as exhibiting "level rigidity" 
0. We will argue here, however, that such descriptions are overly simplistic and apply only in certain energy intervals. 
In fact, since the total number of levels must be conserved, both integrable and chaotic systems exhibit various degrees 
of rigidity depending on the energy scales. 

The main difference between the chaotic and integrable systems lies in the relevant length and energy scales. Chaotic 
systems are characterized by their diffusive behavior Q and the relevant energy scales are the mean level spacing 
A, the Thouless energy E c — D/L 2 , and the inverse scattering time r _1 = D/i 2 , where D <~ vp£ is the diffusion 
coefficient, vf is the Fermi velocity, £ is the mean- free-path, and L is the system size (K = c = 1 below). The Thouless 
scale is associated with the diffusion to the system boundary and can be written as£] E c ~ ^fepA{£/ L), where €f is 
the Fermi energy. 

For chaotic systems, on energy scales E <g; A, the level rigidity statistics A 3 , which is based on the cumulative level 
distribution function (see e.g. M), is 

As ~ § (1) 

indicating the lack of level correlations. The onset of rigidity occurs at E ~ A (and, hence, also known as level 
repulsion) and for A <C E <C E c 

A 3 ~lnf (2) 



This energy range corresponds to length scales between L and y/D/A and the rigidity is related to the diffusive 
character of the multiple traversals of the system. On scales such that E c <C E <C r _1 , which correspond to lengths 
between £ and L and to uncorrelated diffusion in different parts of the system, the level structure is much less rigid 

e. 

E IE 

A 3 „ — ~ -- 3 

E c g A 

where g ~ E c / A 3> 1 is the dimensionless conductance. The existence of the above regimes is, nonetheless, in complete 
agreement with the global rigidity (level number conservation) which can be expressed in the following form: 

djj {8v (e) Sv (e + u)) = (4) 

3 

where v (e) is the level density, 8v (e) = v (e) — (v) and (v) the mean level density. (Here e is the "running" value of 
energy which below will be used interchangeably with ep). This relationship has been discussed, for instance in [|J, 



1 In what follows, we consider only 2D systems. 
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but without taking account of diffusion modes. The diffusion modes are only known to be treated perturbatively, but 
even then the above relationship is satisfied as can be easily seen with the help of eq. (32) in [||. 

In integrable systems (|[) must also be satisfied; otherwise even the notion of the mean level density A would be 
dubious. However, the only relevant scale in this case is 



e ~ v/7fA (5) 

In what follows, we shall derive the expression for the level correlation function (8v (e) 8v (e + uj)) in integrable systems 
that satisfies (Q) and illustrate it on an example of a rectangle with incommensurate sides (a generic integrable system) . 
Applying this expression to evaluation of the A3-statistics, we will find that for E <C e 




(6) 



and for E ^> e 



that is ^-independent. The latter indicates that the spectrum becomes more rigid at large scales, the fact obviously 
related to (Q). The leading terms in (JsJ) and (Q) were a subject of extensive numerical and analytical study in || and 
||. However, the evaluation in the present paper of the level correlation function over the entire energy range allows 
for the evaluation over the entire energy range of A 3 - and E-statistics (and other characterizations of the spectrum 
§) also. 

II. LEVEL CORRELATION FUNCTION 

We approach the derivation of the level-correlation function in two different ways, which will later be shown to be 
equivalent. First, we utilize the scmiclassical formalism wherein the level density is expressed as a sum over periodic 
orbits H , [0 . In such a formalism, the Fourier transform of the level correlation function is given by 

0(t)=x;^(t-T i ) (8) 
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where Aj and Tj are the periodic orbit amplitudes and periods respectively To simplify the calculation, we will 
express all energies in terms of A so that A will be dropped below. It is known that (see (58) in &) 

<f>(t -> oo) = 1/27T (9) 

and it is clear that <f)(t) = for t < T m [ n , where T m i n (~ in 2D) is the shortest periodic orbit. Based on these 
limiting behaviors, we propose a very simple ansatz, namely, 

*(*)= y 2?r ' \ j5 nin (io) 

u ; L ^ ± mill 

Clearly, this ansatz is applicable in any dimension. 

Since the correlation function should be symmetrical with respect to the u) — * — lu transformation, the FT should 
be symmetrical with respect to the t — > — t transformation. The obvious generalization of (g) would be via the 



substitution 5 (t — Tj) — ► [5 (t — Tj) + 6 (t + Tj)] and the ansatz (10) is generalized be means of 

1 



sign (t T m [ n ) sign (t + T min ) 



2 2 

which satisfies <fi (— t) = <f>(t). With such definition, the correlation function becomes 



(11) 



(Sv (e) Sv (e + w)}=8 (lu) - Sm ( "" Tmin) (12) 



Consider now a rectangle whose sides L\ and L2 are such that 



2 



L{jL\ = a (13) 

is irrational. It is also assumed that a < 1 (this assumption is opposite to the assumption a > 1 in ||). We also 
have 

£f = NA = ——7 (14) 
mA 

where A = L\L 2 is the area and N is the mean number of levels below the Fermi level. With these definitions, the 
shortest periodic orbit is the one with length 2Li and 

2L X 2^/2 Q i/4 

T mi „ = — = — — (15) 

The level correlation function for this integrable system is then 

(^(e)^(e + c))=^n- Sin(27r " /£) (16) 

TTUl 

with the definition 

e = V^F/a 1/4 (17) 

(compare this definition with eqs. (7) and (43) in ||). 

The alternative way of derivation is specific to a rectangle (since it is based on the energy spectrum specific to the 
rectangular box) and involves a slight modification of a derivation for a square in Appendix B of . It is shown there 
that the level density can be written as 

x ( \ ST* cos (kLg) 

mi, ?7i2 — — oo y^-^cr J 

The term mi = rri2 = is excluded. (In fact, this term gives the mean level density [§). In @, k = \/2me, kL > 1 
and 



L a = 2a/ (miii) 2 + (m 2 L 2 ) 2 (19) 



Notice that reducing the summation to only positive mi,m2 demonstrates that ( |18| ) is a summation over periodic 
orbits whose length is given by (|l9|). Omitting the rapidly oscillating terms, we find 

(20, 



The sum can be converted to an integral in polar coordinates 

1 



2kL x L 2 



/OO fPmin 
exp (2iAkp) dp —J cos(2Akp)dp 



(21) 



where^j p m i n ~ L\. The appearance of p m m above is related to the fact that the term mi = m2 = is excluded in 
(|l8|). It has the meaning of the length of the shortest periodic orbit, which is in direct relation to ansatz (|Tl]). (Notice 
that when a ~ 1, p m in ~ L\ sa L 2 and the combination of the cut-off and polar coordinates becomes a more tenable 



approximation). Integrating (21) and substituting Afc = Lum/k, we recover the previously obtained result (hq 



2 More consistently, the variable of integration should be a dimensionless L\k <g 1; introducing the latter, however, does not 
affect the final result. 
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III. A 3 - AND E-STATISTICS 



We now proceed to apply ( |l6| ) to the evaluation of A3- and S-statistics. For A 3 -statistics, we begin with 

f°° dt 

A 3 (E)=2j -^(t)G(Et/2) (22) 
G (y) = 1 - F 2 (y) - 3 (F' (y)f , F (y) = siny/y 



(compare with eq. (20) in M). Substituting (11) into (E2J), we find 



A 3 (E) 



E 
2^ 



dt 



G(t) 



BT min /2 



t; 2 



(23) 



This is easily evaluated analytically in terms of algebraic, trigonometric and sine integral functions. However, of 
greater interest are the limiting cases: 



A3 (E < IT* 



e) S — E 
' 15 



1 - 



{Tmi n Ey 
1447T 



and 



A3 (E » T~l ~ e) 



7rT n 



3 (T m i n i?)^ 



For the rectangular, upon substitution of dlq), the first term in (J25|) becomes 



(red) 



(S»e) 



1 



27r 3/2 a l/4 



(24) 



(25) 



(26) 



Remarkably, for a ~ 1, the constant in front of the radical is close (~ 5 percent) to the value in (45) of jg] which 
was obtained by exact summation over all periodic orbits. This supports the utility of the ansatz (Fyl) for the level 
correlation function. 

For completeness we give the results for S-statistics as well. In the same limiting cases, we find 



S (E « ~ e) 



E 



1 



(T m i n E) 



and 



S (£ » TTi ~ e) s 



7rT n 



sin(r min S) 2cos(T min £;) 



mm 



(Tmin-E) 



(27) 



(28) 



We observe that the accuracy of A 3 in the limiting cases is higher than that of S. This is because A 3 describes the 
cumulative behavior of the levels in the spectrum and is thus a more appropriate characteristic for the rigidity. We 
also observe that the last two terms in (j^) are rapidly oscillating. The most interesting feature of both statistics for 
large energy intervals is that the leading terms are ^-independent (but do depend on e). 



IV. DISCUSSION 

The central result of this work is (|l2|) . This simple form of the correlation function is in excellent agreement with the 
limiting behavior of the A3-statistics for both small and large energy scales and is consistent with the overall spectral 
rigidity (|]). Having a closed- form correlation function, applicable at all energy scales, allows for the evaluation of the 
closed-form expressions for A 3 - and S-statistics as well, with the corrections in the appropriate regimes given by eqs. 
@)-@ and eqs. @-(||) respectively. 



The significance of the corrections to the leading terms in (24) and ( |27| ) is that they describe the onset of level 
correlations at small scales, which are dominated by the 5 -function in ([l2]) (corresponding to the absence of correlations 
in the zeroth order). The corrections to the leading terms at large scales, (p5|) and (p8|), are also of interest especially 
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in view of the issues related to the symmetry-breaking perturbations (see below). For verification purposes, we are 
presently working on the numerical evaluation of S-statistics which will be reported elsewhere. 

We emphasize that (12) was obtained both via ansatz (]ll]) and, for the specific case of a rectangle with incommen- 
surate sides, by a direct evaluation in the energy space. In both circumstances the small inaccuracy is related to the 
continuous nature of the approximation used versus the discrete behavior in the limit of short orbits. 

Clearly, if a term 56 (t), such that 56 (T < T^ in ) = and 56 (±oo) = 0, is added to (|lTJ) , eq. (Q) will still be satisfied 
and the leading terms in (24) and (pTj), determined by the (5-function in (p"2|), will remain unchanged. Furthermore, 
since the relevant scales are TL, n > T m i n , the corrections introduced at large energy scales will have the same functional 
form as, and will be of the order of, ( p5| ) and (28). The question of 5<p(t) may be relevant in discussion of the effect 
of time-reversal symmetry-breaking terms, such as magnetic field. We hope to address this question in a future 
publication. 
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